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By applying density functional theory (DFT) approximations, we present a first- 
principles investigation of elastic properties for the experimentally verihed phases of 
a metallic perovskite LaNiOa. In order to improve the accuracy of calculations, at 
hrst we select the most appropriate DFT approaches according to their performance 
in reproducing the low-temperature crystalline structure and the electronic density 
of states observed for the bulk LaNiOa. Then, we continue with the single-crystal 
elastic constants and mechanical stability for the most common rhombohedral as 
well as high-temperature cubic and strain-induced monoclinic phases. Together with 
the calculated single-crystal elastic constants, the deduced polycrystalline properties, 
including bulk, shear, and Young’s moduli, Poisson’s ratio, Vickers hardness, sound 
velocities, Debye temperature, and anisotropy indexes, remedy the existing gap of 
knowledge about the elastic and mechanical behaviour of LaNiOa, at least from a 
theoretical standpoint. 
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I. INTRODUCTION 


Among the remarkable family of perovskite oxides, lanthanum nickelate LaNiOa is a rare 
example characterized by paramagnetic metallic behaviour down to the lowest tempera¬ 
tures and being structurally compatible with many active functional layers [1, 2], For a 
long time, the technological interest in this material has been limited to the development 
of highly conductive bottom electrodes for various ferroelectric thin-hlm devices [3-5]. But 
recently, the enthusiasm was renewed by experimental recognition of electric-held tunable 
metal-to-insulator transition and theoretical hndings that heterostructures and superlattices 
composed of thin LaNiOs layers separated by wide-gap insulators could possibly exhibit 
high-temperature superconductivity [6-8]. Of course, the fact that the hunt for novel super¬ 
conducting materials is at the vanguard of condensed matter physics research sheds a new 
light on the outlook of LaNiOs, but on the other hand, it does not automatically mean that 
the bulk properties of this oxide are already fully explored. For example, although one can 
hnd a variety of theoretical studies [2, 9-12] based on conventional density functional theory 
(DFT) and beyond-DFT techniques in which the electronic structure of LaNiOs is exten¬ 
sively investigated and compared to the experimental data, in the literature scarcely any 
hrst-principles calculations on the elastic properties of LaNiOs have been reported so far. 
The situation is not much better with the experiment, especially regarding the single-crystal 
measurements. Having in mind that the list of potential applications of LaNiOs continues 
to grow [13-15], it becomes obvious that the deep knowledge of elastic features could be 
useful for both technological and fundamental reasons. 

The main goal of our work is to evaluate single-crystal elastic constants, polycrystalline 
moduli, and other mechanical properties of LaNiOs using first-principles calculations. In 
order to select the most appropriate DFT-based approaches for it, at hrst we focus attention 
on the geometrical parameters and electronic structure of rhombohedrally distorted phase 
(space group R3c, No. 167) which is observed in a wide range of temperatures. Then, we 
continue with the elastic parameters and extend our study to the high-temperature cubic 
phase {Pm3m, No. 221) [16] together with monoclinic {C2/c, No. 15) [17] structure that is 
available under tensile or compressive strains induced by various substrates. We hope that 
the obtained data will be helpful for many practical applications related to the mechanical 
behaviour of bulk LaNiOs, thin hlms, and heterostructures. 
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II. COMPUTATIONAL DETAILS 


All the calculations reported in this paper were performed with the CRYSTAL14 program 
[18] which expands the crystalline orbitals as linear combinations of atom-centred Gaussian- 
type functions. A small-core Hay-Wadt pseudopotential [19] was used for describing the core 
electrons (ls^2s^2p®3s^3p®3(i^°4s^4p®4(i^°) of La atom, while the valence part (5s^5p®5(i^6s^) 
of this basis set was adopted from LaMnOa study [20]. Concerning the description of Ni and 
O atoms, the all-electron triple-zeta plus polarization quality basis sets, specihcally derived 
for periodic solid-state treatment, were taken from [21]. 

For the sake of simplicity, we have modelled the paramagnetic behaviour of LaNiOs by 
non-magnetic calculations. The default values of parameters that dehne the convergence 
threshold on total energy and truncation criteria for bielectronic integrals were set to more 
severe ones: (8) and (8 8 8 8 16), respectively. In addition, while performing full geometry 
optimization the allowed root-mean-square values of energy gradients and nuclear displace¬ 
ments were correspondingly tightened to 0.00006 and 0.00012 a.u. In order to improve the 
self-consistence held convergence, the Kohn-Sham matrix mixing technique (at 80%) to¬ 
gether with Anderson’s method [22] were applied. The shrinking factor of 32, 16, and 8 
was used for reciprocal space sampling, resulting in 969, 417, and 150 independent k points 
in the hrst irreducible Brillouin zone for cubic, rhombohedral, and monoclinic phases of 
LaNiOs, respectively. All the other important settings for calculations were left at their 
default values, which in turn can be found in CRYSTAL14 user’s manual [23]. 

A variety of DFT approaches was tested in order to hnd out which of them exhibit 
the closest resemblance to the low-temperature experimental data on crystalline structure: 
starting from standard local density (LDA [24, 25]) and generalized gradient (PBE [26]) 
approximations up to revised functionals for solids (PBEsol [27], SOGGA [28], and WG 
[29]) and hybrids (mBlWG, PBEO [30], and HSE06 [31]). Note that mBlWG stands for 
modihed BIWG approach [32] since instead of PW functional [33] we have employed the 
correlation part from PBE framework. The revised functionals, SOGGA and WG, were also 
combined with PBE correlation. 

A fully automated procedure ELASTGON [34, 35] developed for calculating single-crystal 
elastic constants Cij has been applied for all three phases of LaNiOs. For each independent 
strain, two deformed structures with the magnitude of the strain reaching -0.01 and 0.01 
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were considered. The internal atomic positions were allowed to relax with the deformed 
cell shape and volume remaining hxed. Polycrystalline elastic properties, starting from the 
bulk and shear moduli, were derived from the computed single-crystal elastic constants Cij 
using the Voigt [36] and Reuss [37] approximations. The mean values of these parameters 
calculated via the averaging method proposed by Hill [38] have been utilized to evaluate 
Young’s modulus, Poisson’s ratio, sound velocities, and Debye temperature according to 
well-known expressions that can be found elsewhere [39, 40]. For comparison purposes, the 
bulk modulus, denoted as Rbm, has been obtained by htting the total energy as a function 
of volume to the third-order Birch-Murnaghan (BM) equation of state [41]. The htting was 
done at 11 points within the range of 0.92-1.08 variation of the initial volume by another 
automated procedure, named EOS, which is also implemented in the CRYSTAL14 code. 


III. RESULTS AND DISCUSSION 
A. Performance of DFT approaches 


Since many specihc material properties like elastic constants, sound velocities, and Debye 
temperature are critically dependent on change in geometry of the system, we primarily 
focus our attention on the accuracy of the calculated lattice constants and other structural 
parameters that are shown in Fig. 1. The obtained values together with low-temperature 
experimental data for rhombohedral LaNiOa are presented in Table I. The mean absolute 
relative error (MARE) was evaluated according to the expression 
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in which and are the calculated and experimental values of the considered param¬ 
eter, respectively. Note that we did not apply the zero-point anharmonic expansion (ZPAE) 
corrections for the experimental data. On one hand, the ZPAE corrections are straightfor¬ 
wardly applicable only for the cubic systems [42]. On the other hand, ZPAE can expand 
the equilibrium lattice constant by 1% for light atoms like Li and much less for heavy atoms 
[43], thus it should not have a noticeable influence on material like LaNiOs. What is more, 
we did not analyse the geometry of cubic and monoclinic phases of LaNiOs because of the 
thermal expansion of the former and internal strains of the latter meaning that the direct 
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FIG. 1. The crystalline structure of rhombohedral {R3c) LaNiOs. The structural parameters are 
presented in Table I. 

comparison between experiment and fully-relaxed zero-temperature calculations would not 
be accurate. 

As it could be expected, Table I reveals that LDA and PBE functionals tend to un¬ 
derestimate and overestimate lattice constant a together with angle a, respectively. As a 
direct consequence, this tendency is also reflected in a deviation of the volume V values. 
However, the following results are in a good agreement with the recent plane-wave calcu¬ 
lations [2, 12] demonstrating the reliability of our chosen basis set. Although the revised 
functionals - PBEsol, SOGGA, and WG - appear worse than PBE for the evaluation of 
the Ni-O-Ni angle, their overall performance is systematically improved, especially for the 
aforementioned structural parameters, a, a, and V. The decrease of MARE from 1.15% for 
PBE to 0.74% for SOGGA and impressive numbers of 0.41 and 0.35% for PBEsol and WG, 
correspondingly, leaves no doubt about the effectiveness of the modihcations made to the 
revised density functionals for solids. A slightly worse result of SOGGA could be influenced 
by its different analytical form or more tightly bounded exchange [28]. 

One can note that the inclusion of the standard portion (25%) of the exact Hartree- 
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TABLE I. Calculated structural parameters of rhombohedral LaNiOs compared to the experimental 
data [54] at 1.5 K. Lattice constant a and bond distance Ni-0 are given in A, volume V is given 
in A^, angles a, Ni-O-Ni, and O-Ni-0 are given in degrees. MARE (in %) stands for the mean 
absolute relative error. The numbers in brackets (in %) represent absolute relative errors for each 
of the considered parameters. 



LDA 

PBE 

PBEsol 

SOGGA 

WG 

inBlWG 

PBEO 

HSE06 

Expt. 

a 

5.332 

5.444 

5.375 

5.361 

5.379 

5.363 

5.396 

5.399 

5.384 


(0.95) 

(1.11) 

(0.16) 

(0.42) 

(0.09) 

(0.39) 

(0.23) 

(0.27) 


a 

60.68 

61.09 

60.89 

60.83 

60.90 

60.42 

60.44 

60.49 

60.86 


(0.30) 

(0.37) 

(0.05) 

(0.05) 

(0.06) 

(0.73) 

(0.69) 

(0.60) 


V 

108.86 

116.85 

112.03 

110.99 

112.25 

110.07 

112.20 

112.49 

112.48 


(3.22) 

(3.89) 

(0.40) 

(1.32) 

(0.20) 

(2.14) 

(0.25) 

(0.01) 


Ni-0 

1.904 

1.958 

1.926 

1.919 

1.927 

1.908 

1.922 

1.924 

1.933 


(1.50) 

(1.29) 

(0.37) 

(0.73) 

(0.31) 

(1.27) 

(0.58) 

(0.45) 


Ni-O-Ni 

169.05 

164.77 

167.05 

167.60 

167.04 

170.70 

169.92 

169.38 

164.82 


(2.57) 

(0.03) 

(1.35) 

(1.69) 

(1.35) 

(3.57) 

(3.09) 

(2.77) 


O-Ni-0 

89.16 

88.58 

88.88 

88.96 

88.88 

89.46 

89.40 

89.33 

88.78 


(0.43) 

(0.22) 

(0.12) 

(0.21) 

(0.12) 

(0.77) 

(0.71) 

(0.63) 


MARE 

1.49 

1.15 

0.41 

0.74 

0.35 

1.48 

0.92 

0.79 



Fock (HF) exchange - as it is in the PBEO scheme - has a positive impact on the overall 
performance of PBE functional, since MARE now gets reduced to 0.92%. Eurthermore, 
the rejection of the long-range HF-type exchange in the HSE06 framework seems to have 
even a better effect. As HSE06 outperforms PBEO for nearly every structural parameter of 
rhombohedrally distorted LaNiOs, its MARE takes a value of 0.79% which is very close to 
the result of SOGGA. Interestingly, HSE06 functional manages to almost ideally reproduce 
the experimentally determined value of volume V, although neither lattice constant a nor 
angle a alone are improved in comparison to what was obtained using PBEsol or WG. In 
this case, HSE06 benehts from an error cancellation which, according to the expression 

V = o^a/I -f- cos^ q;( 2 cosa — 3), (2) 
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is likely if a and a are correspondingly overestimated and underestimated (or vice versa) in 
respective proportions. 

Another popular hybrid functional for the solid-state calculations, named BIWC, incor¬ 
porates 16% of the exact full-range HF exchange. In this work, a modihed version mBlWC 
with PBE correlation was applied in order to find out the effect of HF exchange on solely the 
WC framework. Unfortunately, Table I apparently indicates that mBlWC not only degrades 
the accuracy of its predecessor WC for all the structural parameters, but also is the worst 
performer among the considered hybrids with MARE reaching 1.48%. It becomes obvious 
that the inclusion of some part of HF exchange is harmful for WC and, most probably, 
PBEsol and SOGGA functionals which already demonstrate reliable results in reproducing 
the crystalline structure of rhombohedral LaNiOs. 

On the whole, the evaluated mean discrepancies between different DFT calculations and 
experiment do not exceed 1.5%, and in most cases this limit may seem to be acceptable. 
However, for the sake of greater accuracy, we tighten functional selection criterion to 1%, 
therefore LDA, PBE, and mBlWC appear to be the ones that should be eliminated before 
making further calculations of elastic properties. 

But this is not the end of the story. Since simply reproducing the correct ground-state 
crystalline structure is not a sufficient criterion to properly evaluate the quality of the DFT 
approach performance [2], in Fig. 2 we also compare the density of states (DOS) obtained 
using the remaining functionals with experimentally observed valence band of LaNiOs. Here, 
it is clearly seen that hybrids fail to correctly describe the electronic structure exhibiting 
only two signihcantly shifted spectral features, while revised functionals demonstrate a fairly 
good agreement with the main energy peak positions A, B, and C. Needless to say, such 
combination of accurately reproduced electronic and crystalline structures shown by PBEsol, 
SOGGA, and WG automatically captures our sympathy. Besides, if one would take into 
account the computational time and resources required for the calculations, they would 
have another advantage over the hybrid scheme of PBEO and HSE06. Thus, to sum up, it 
can be stated that the revised functionals for solids appear to be the best choice among the 
considered DFT approximations for the next calculations of elastic properties of LaNiOs. 
And we hope that the quality of the description of the other two phases will remain at the 
same level as it is now for rhombohedral LaNiOs. 
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FIG. 2. The comparison of total DOS calculated using (a) hybrid and (b) pure DFT approaches 
with (c) X-ray photoelectron spectroscopy LaNiOa spectrum [11]. The Fermi energy is set at zero. 


B. Elastic properties 

Single-crystal elastic constants Cij of LaNiOs calculated by applying revised functionals 
are presented in Table II. Unfortunately, no experimental or theoretical data are currently 
available to deal with, therefore we can only make a comparison between our own results. 
From Table II it is seen that PBEsol, SOGGA, and WG perform quite similarly, with SOGGA 
on average giving slightly higher values than PBEsol and WG. This indicates that the ob¬ 
tained results are reliable and could be useful for various practical applications. Goncerning 
the elastic properties of LaNiOs, one can note that rhombohedral phase is more resistive 
against uniaxial deformation along the z direction, while monoclinic phase is harder to de¬ 
form along the x axis. The average value of all three functionals for and reaches 
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TABLE II. Single-crystal elastic constants Cij (in GPa) calculated for different phases of LaNiOs. 


Cubic (PmSm) Rhombohedral (i?3c) Monoclinic (C^ljc) 



PBEsol 

SOGGA 

WC 

PBEsol 

SOGGA 

WG 

PBEsol 

SOGGA 

WC 

Gii 

405.5 

415.2 

405.3 

264.6 

268.6 

266.2 

341.9 

348.1 

342.9 

C12 

119.2 

120.0 

120.7 

161.6 

165.1 

162.4 

156.3 

164.6 

159.3 

Gi3 




150.1 

151.3 

151.2 

136.5 

139.6 

137.3 

Gi5 




17.6 

17.6 

16.9 

-3.9 

-3.7 

-3.8 

C22 







269.3 

270.4 

270.7 

C23 







145.1 

148.3 

148.6 

C25 







16.7 

17.2 

16.9 

G 33 




341.1 

348.9 

342.2 

301.1 

307.2 

302.8 

G 35 







-30.2 

-30.9 

-29.9 

C44 

103.6 

105.5 

103.5 

98.8 

101.6 

99.5 

46.2 

45.9 

46.9 

C46 







0.7 

-1.7 

- 1.2 

C55 







82.1 

83.9 

82.8 

Gee 




51.5 

51.7 

51.9 

104.9 

107.8 

105.7 


~ 344 GPa. The resistance against uniaxial compression or tension along the y direction 
is found to be practically the same for both structures. Although cubic phase exhibits the 
highest incompressibility of all under the uniaxial strains with = € 22 '°' = ^' 33 ^'^' ~ 409 
GPa, its resistance against shear deformations is reduced to ~ 104 

GPa. A similar value is demonstrated by rhombohedral structure for shear deformations 
around the x and y axes ~ 100 GPa) and by monoclinic phase for shear 

deformation around the 2 : axis ~ 106 GPa). But despite that, the latter is much 

easier to deform around the x axis with ~ 46 GPa, whereas the former - around 

the 2 : axis with ^ 52 GPa. Some of the elastic constants, such as or 

are negative and this may imply that monoclinic phase of LaNi 03 could be mechanically 
unstable. However, this is not the case, since the stability restrictions on Gjj, following from 
the requirement that the strain energy must always be positive for all applied strains, are 
fully satisfied. It indicates that cubic, rhombohedral, and monoclinic structures of LaNi 03 
are mechanically stable at zero temperature and zero pressure conditions. The appropriate 
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stability criteria on Cij for different phases of material can be found elsewhere [44, 45]. 

The knowledge of single-crystal elastic constants Cij allows one to evaluate mechanical 
and other important properties of polycrystals. The list of the parameters which describe the 
polycrystalline behaviour of LaNiOs is given in Table III. For the sake of convenience, Table 
III presents values of bulk and shear moduli obtained by applying not only the averaging 
method recommended by Hill (H, G), but also the initial isostrain and isostress approxima¬ 
tions correspondingly suggested by Voigt (Hy, Gy) and Reuss (Hr, Gr). It is worth to note 
a very good agreement between the bulk modulus B and the one evaluated from the BM 
equation of state Rrm indicating a high numerical accuracy of our calculations. Regarding 
the available experimental data, a value of ~ 1^8 GPa from the room-temperature 

measurements [46] nicely corresponds to ^ 200 GPa, especially if one would take 

into consideration thermal effects which tend to reduce the elastic moduli. A comparison 
between the different phases of LaNiOs in Table III reveals that all three structures exhibit 
a similar level of compressibility because the difference in their bulk moduli B does not 
exceed ~ 8 %. However, the cubic structure is more resistive to shear and uniaxial deforma¬ 
tions: 119 GPa against ~ 73 GPa and G'^°°' ~ 74 GPa, and ^ 301 

GPa against 194 GPa and si 197 GPa. This fact is also reflected in the 

lower value of Poisson’s ratio and higher values of sound velocities and Debye temperature 
compared to the rhombohedral and monoclinic phases. The available experimental data for 
Debye temperature = 420 K [47] indicate a fairly good agreement with calculations 

^Rhomb. K, since the deviation from experiment does not exceed ~ 13%. Further anal¬ 

ysis of elastic parameters shows that the Vickers hardness for cubic LaNiOs is more than two 
times larger in comparison to the rhombohedral and monoclinic LaNiOs: fa 13.7 GPa 

against 6.1 GPa and 6.2 GPa, respectively. This well-known mechanical 

property, dehned as the ability of a material to resist against being dented or scratched 
by another, was obtained according to the expression By = 0.92(G/R)^'^^’^G°™® recently 
proposed by Tian et al. [48]. Interestingly, the apparent distinction in hardness does not 
change the ductile nature of cubic, rhombohedral, and monoclinic LaNiOs, although the 
latter two structures demonstrate a more pronounced ductility with corresponding B/G 
ratios of ~ 2.75 and ~ 2.71 compared to ~ 1.82 of the cubic phase. The most common 
critical value that separates ductile and brittle character was introduced by Pugh [49]: if 
B/G > 1.75, the material behaves in a ductile manner, otherwise, its behaviour should be 


10 



TABLE III. Polycrystalline bulk modulus (in GPa) from Voigt (By), Reuss (Br), and Hill (B) 
approximations, bulk modulus Rbm (in GPa) and its pressure derivative from BM equation of 
state, shear modulus (in GPa) from Voigt (Gy), Reuss (Gr), and Hill (G) approximations. Young’s 
modulus E (in GPa), Poisson’s ratio longitudinal sound velocity ul (in m/s), transverse sound 
velocity ut (in m/s), average sound velocity v (in m/s), Debye temperature 0 d (in K), Vickers 
hardness Hy (in GPa), B/G ratio, anisotropy index in compressibility As (in %), anisotropy 
index in shear Aq (in %), and universal anisotropy index A^ calculated for different phases of 
LaNiOs. 


Cubic {Pmim) Rhombohedral (i?3c) Monoclinic (C^ljc) 



PBEsol 

SOGGA 

WC 

PBEsol 

SOGGA 

WC 

PBEsol 

SOGGA 

WC 

By 




199.3 

202.4 

200.5 

198.7 

203.4 

200.8 

Rr 




197.5 

200.6 

198.6 

196.7 

201.2 

199.0 

B 

214.7 

218.4 

215.6 

198.4 

201.5 

199.6 

197.7 

202.3 

199.9 

Rbm 

214.7 

218.5 

215.6 

198.9 

202.5 

200.0 

199.3 

201.9 

199.3 

Rbm 

4.26 

4.28 

4.25 

4.28 

4.26 

4.28 

4.22 

4.26 

4.22 

Gy 

119.4 

122.3 

119.0 

77.1 

78.8 

77.5 

78.3 

79.1 

78.5 

Gr 

116.5 

119.0 

116.2 

66.7 

68.0 

67.6 

69.1 

69.1 

69.2 

G 

118.0 

120.7 

117.6 

71.9 

73.4 

72.6 

73.7 

74.1 

73.9 

E 

299.1 

305.7 

298.5 

192.6 

196.5 

194.2 

196.6 

198.0 

197.3 

V 

0.27 

0.27 

0.27 

0.34 

0.34 

0.34 

0.33 

0.34 

0.34 

WL 

7128 

7168 

7140 

6368 

6393 

6396 

6387 

6411 

6419 

Vt 

4014 

4043 

4012 

3148 

3166 

3165 

3187 

3180 

3194 

V 

4466 

4497 

4465 

3534 

3554 

3553 

3575 

3569 

3584 

Su 

595.9 

601.8 

595.3 

470.3 

474.3 

472.4 

475.7 

476.3 

476.5 

Hy 

13.65 

13.96 

13.50 

5.99 

6.11 

6.05 

6.29 

6.19 

6.24 

B/G 

1.82 

1.81 

1.83 

2.76 

2.75 

2.75 

2.68 

2.73 

2.71 

Ab 

0.0 

0.0 

0.0 

0.45 

0.45 

0.48 

0.51 

0.54 

0.45 

Ag 

1.23 

1.37 

1.19 

7.23 

7.36 

6.82 

6.24 

6.74 

6.30 

RU 

0.12 

0.14 

0.12 

0.79 

0.80 

0.74 

0.68 

0.73 

0.68 
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related with brittleness. 

Elastic anisotropy can have considerable effects on durability of the material through the 
formation of microcracks, therefore this parameter plays an important practical role not less 
than the others. The percentage anisotropy in compressibility and shear for polycrystals, 
originally suggested by Chung and Buessem [50], is dehned as Ab = l00{By—Bji)/{B y+Bb) 
and Ag = 100(Gv — Gr)/(Gv + Gr); respectively. Here, a value of zero represents elastic 
isotropy, while a value of 100% identihes the largest possible anisotropy. From Table III it is 
seen that rhombohedral and monoclinic phases of LaNiOs exhibit a similarly weak anisotropy 
in compressibility with yl^homb. ^ 0.46% and Ri 0.5%, but their anisotropy in shear 

appears to be more pronounced, as the numbers accordingly increase up to Ri 7.14% 

and r^ 6.43%. In the meantime, the perfect isotropy in compressibility for the cubic 

phase accompanied by its relatively weak anisotropy in shear A^^' Ri 1.26% speaks for the 
highest degree of elastic isotropy among the considered structures of LaNiOs. And indeed, 
the following Ending is confirmed by the universal anisotropy index = 5(Gv/Gr) — 
{By/Bb) — 6 [51], the values of which can also be found in Table III. By taking into account 
both the bulk and the shear contributions, A^ represents a universal measure to quantify 
the single-crystal elastic anisotropy. For an isotropic crystal, Al^ = 0 and deviations from 
zero determine the extent of elastic anisotropy. Gladly, the comparison of calculated values 
^Cub. ~ 0.13, Ri 0.78, and Amoii. ~ 0.7 shows a full compatibility with our previous 

observations drawn from As and Aq- 

In general, one can notice that polycrystalline parameters derived for rhombohedral and 
monoclinic phases of LaNiOs bear a very close resemblance to each other, since in most 
cases the discrepancy does not exceed ~ 10%. On one hand, it implies that transformation 
of rhombohedral LaNiOs into monoclinically distorted structure induced by compressive or 
tensile strains does not have a strong influence on its mechanical behaviour. But on the 
other hand, one should take into consideration that here we are dealing with a fully relaxed 
monoclinic system which is slightly different from non-relaxed ones experimentally identi¬ 
fied on SrTiOs and LaAlOs substrates [17]. However, we believe that our elastic parameters 
pretty reliably describe monoclinic LaNiOs as long as the thickness of the film is sufficiently 
large to prevent from the occurrence of metal-to-insulator transition. Experimental investi¬ 
gation of electrical resistivity [52] and electronic structure [53] has revealed that films that 
are thick enough (10 nm or somewhat less depending on the substrate) show clearly pro- 
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nounced metallic behaviour that is typical for bulk LaNiOs and can be properly represented 
at DFT level. But for the ultrathin LaNiOs dims composed of only several monolayers, a 
metal-to-insulator transition is observed and a standard DFT scheme fails to reproduce it 
[53]. Thus, in order to accurately describe the elastic properties of ultrathin LaNiOs struc¬ 
tures, one should apply more sophisticated theoretical methods. In other cases, our given 
data seem to be sufficient. 


IV. CONCLUSIONS 

Elastic properties allow one to evaluate many important mechanical and thermodynamic 
features of the material, thus in this paper we remedy the existing gap of knowledge by 
presenting single-crystal elastic constants and deduced polycrystalline parameters for the 
metallic perovskite LaNiOs. For the sake of greater accuracy, a bunch of DFT approaches 
has been tested to reproduce the experimentally observed low-temperature crystalline and 
electronic structures of bulk LaNiOs. As the revised functionals for solids have demon¬ 
strated the best performance, they were chosen for the next calculations of elastic properties 
for the most common rhombohedral as well as high-temperature cubic and strain-induced 
monoclinic phases of LaNiOs. The obtained results indicate that all three structures are 
mechanically stable, behave in a ductile manner, and exhibit a similar level of compressibil¬ 
ity. However, rhombohedral and monoclinic LaNiOs are softer against uniaxial and shear 
deformations in comparison to the cubic phase, and this fact is somewhat reflected in the 
higher values of the Poisson’s ratio and lower ones for the Vickers hardness, sound velocities, 
and Debye temperature. Concerning the elastic isotropy, it can be stated that, in general, 
LaNiOs is more anisotropic in shear than in compression but the overall anisotropy is not 
strongly pronounced. It is also worth to note that rhombohedral and monoclinic phases 
demonstrate a very close polycrystalline behaviour which in turn implies that LaNiOs re¬ 
tains its mechanical properties when transforming from thin him to bulk material. We expect 
the obtained data and proposed insights to be useful for various practical applications and 
possible future experiments. 
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